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The influence of a constant homogeneous external magnetic field H on the formation and stability 
of quark droplets is investigated within a simple Nambu - Jona-Lasinio model by using a thermo- 
dynamic approach. For a vanishing magnetic field stable quark droplets, which are schematically 
the bags of massless quarks, are allowed to exist only at G > Gbag, where G is the quark coupling 
constant, Gb ag = l-37G cr u, and G cr it is the value of the coupling constant above which chiral sym- 
metry is spontaneously broken down. On the other hand, a nonvanishing external magnetic field 
£f*l | can induce the stability of quark droplets so that they may exist even at G < Gbag- In this case, 

depending on the value of H, quark droplets are composed either of massive or massless quarks. 
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I. INTRODUCTION 

> ' 

' At present time, it is generally accepted that QCD is the theory of strong interactions between mesons and baryons. 
Moreover, it should describe the physical phenomena in dense hadronic matter as well. Since the usual (perturbative) 
t— I ■ methods of QCD are applicable only at sufficiently high values of the baryon density p, in the region of low and 
moderate densities one should use nonperturbative methods, e.g. lattice QCD calculations 1], or effective field theory 
considerations. The latter ones are especially useful if the influence of external magnetic field on the properties of 
dense hadronic matter is investigated (in this case the incorporation of strong magnetic field into lattice calculations is 
not eleborated sufficiently). It is well-known that at low (baryon) densities p < po, where po = 0.16 fm~ 3 is the normal 
density of baryons inside nuclei, hadronic matter is just the ordinary nuclear matter. The widely-known relativistic 
effective theory for nuclei and nuclear matter is the Walecka model 0, • In this model the protons and neutrons, 
i.e. the simplest baryons (nucleons), are point particles. Moreover, the existence of stable objects with large baryonic 
, number, which are droplets of simplest massive baryons and hence might be identified with nuclei, is predicted in the 
1 framework of it with the help of the condensed matter (thermodynamic) approach. 

Unfortunately, the Walecka model is unrealistic in the region of moderate densities near to the chiral phase transition, 
where it is necessary to take into account the internal quark structure of protons and neutrons. Similar conditions 
exist in the experiments on relativistic ion-ion collisions, where it is expected that for a short time stable clumps of 
quark matter (usually called quark droplets) with rather large baryonic number might be created. Another example 
is presented by neutron stars, where under extreme pressure nucleons might fuse together, forming droplets of quarks 
(evidently, the baryonic density inside quark droplets is larger, than p ). 

One of the widely-known effective theory of quark- hadronic matter is the Nambu - Jona-Lasinio (NJL) model 4] , 
which was successfully used in describing low-energy meson physics |5| , diquarks [6J, [3 and meson-baryon interactions 
0,0. In the framework of the NJL model quarks are treated as elementary point-like particles, whereas mesons are 
collective quark-antiquark exitations of its ground state. This model can also serve as a starting point for investigations 
of the mass spectrum of the octet and decuplet of baryons (see, e.g. pLOL IllL IT3 | . However, in order to be applicable to 
dense quark-hadronic matter, the NJL model should predict the existence of stable multiquark droplets as well, i.e. of 
stable finite-size objects in space, which are composed of more than 3 quarks and do not tend to shrink to a point or to 
expand over the whole space. One particular way to investigate such objects in the NJL model and which we will follow 
in this paper is based on the use of the general condensed matter approach restricted to vanishing temperature T. In 
this case a stable multiquark object of the system is identified with a droplet of the dense quark phase, surrounded 
by the vacuum phase with zero baryonic density. Hence, if there is the thermally equilibrated coexistence of these 
two phases, then one may conclude about the existence of stable quark droplets. Just in this way, the problem of 
quark matter stability was investigated at zero external magnetic field H in |L3L Il4l Il5| in the framework of a chirally 
symmetric NJL model. Let us quote the two main results of these papers: i) at H = 0, stabilization of (multi)quark 
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droplets is only possible, when they consist of massless quarks, ii) the stability requirement puts strong constraints 
on the acceptable values of model parameters (coupling constant G, cutoff parameter A etc). In particular, it was 
shown that in the thermodynamic approach based on the NJL model, quark droplets can be schematically interpreted 
as MIT-type bags. However, in the MIT 3-quark bag model the bag constant is a phenomenological input parameter 
|16| , whereas in the NJL model it is a dynamical quantity originating from the dynamical breaking of chiral symmetry. 

Taking into account the fact that inside neutron stars the magnetic field strength might reach extremely high 
values (the influence of ultra-strong magnetic fields on physical processes inside neutron stars and on properties of 
dense quark matter are under permanent investigations |171 Il8t Il9l |20 | h it seems to us important to reinvestigate 
the problem of the quark droplet stability in the framework of NJL-type models at H ^ 0. In the present paper, we 
study, in particular, the influence of a constant homogeneous external magnetic field H in a simple NJL model with 
two quark flavors. In Section 2, the stability constraints on the model parameters as well as on the chemical potential 
/i are presented. Apart from this, the phase structure of the model is found at T = 0, H = and \i ^ 0. On this 
basis, it is shown that at H = stability of multiquark droplets, which consist here only of massless quarks, occurs 
at G > Gbag = l-37G cr it with G cr it = it 2 /(3NfA 2 ) and Nf being the number of flavored quarks. In Sections 3, 4 and 
5 the problem of the existence of multiquark droplets is discussed at H ^= 0. In particular, it is argued on the basis 
of our earlier papers |2ll |22| that for each value G smaller than Gbag, the stability of quark droplets is induced by 
an external magnetic field. Moreover, for each G £ (G cr it, Gbag) it is shown that in contrast to the case H = 131, 
there exists an external magnetic field value H a ~ 10 19 Gauss such, that at < H < H a multiquark droplets may be 
formed from massive quarks. Moreover, we have shown using energetical arguments that the magnetic field promotes 
the creation of stable quark droplets and enhances the stability of it GVGI1 clt G > Gb a g ■ 



II. QUARK DROPLET STABILIZATION AT H = 
A. NJL model and stability constraint 



Let us first give several (very approximative) arguments motivating the chosen structure of our QCD-motivated 
NJL model introduced below. For this aim, consider two-flavor QCD, Nf = 2, with nonzero chemical potential 
and color group SU C (N C ) (in our case N c = 3). By integrating in the generating functional of QCD over gluons 
and further "approximating" the nonperturbative gluon propagator by a 6— function, one arrives at an effective local 
chiral four-quark interaction of the NJL type describing low energy hadron physics. Finally, by performing a Fierz 
transformation of the interaction term and taking into account only scalar and pseudo-scalar (qq) interactions, one 
obtains a four-fermionic model given by the following Lagrangian (in Minkovski space-time notation) 1 

L = q[iYd» + VJ°]q + G[(qq) 2 + (qij 5 rq) 2 }. (1) 

In Q /i > is the chemical potential of quarks which is the same for all flavors. Furthermore, we use the notations 
t = (t 1 , r 2 , t 3 ) for Pauli matrices in the flavor space. Clearly, the model is invariant under the continuous chiral 
SU(2)l x SU(2)r group. Apart from this, it is invariant under the discrete chiral transformation, q —>■ i^ b q. (Since 
the baryonic charge is conserved, there exists also a corresponding t/(l)v-invariance of the model). At tree level, the 
Lagrangian \Q contains besides of the chemical potential only one free model parameter, i.e. the coupling constant 
G. Clearly, when including quantum effects (loops) and performing necessary regularizations, an additional cutoff 
parameter A will appear. 

An illustrative way for the investigation of stable quark droplets in quantum field theory is based on a schematic 
analogy with the thermodynamical approach to the liquid-gas first order phase transition. In such a picture one can 
consider the quark droplet as a droplet of the dense liquid phase of the system, whereas the low-density gas phase 
corresponds to the vacuum (the zero density state) which surrounds the quark droplet (strictly speaking, such an 
analogy makes really sense only for systems at finite temperature). Moreover, we suppose that the two phases are in 
the thermal and chemical equilibrium. Of course, in this way, i.e. when ignoring surface effects of droplets, we have 
an adequate consideration of the quark matter stability problem only for quark droplets which contain a very large 
number of quarks inside a large volume (for comparison see j2^|, where surface effects of droplets were taken into 
account). It turns out that not for all values of input parameters of the system this schematic picture of the existence 
of stable quark droplets is realized. Thus, the stability problem of quark droplets arises, which might be expressed in 
terms of several quantitative constraints on the parameters of the system. To solve these constraints means to find 
the values of the parameters of the system at which there might appear stable quark droplets at moderate densities. 
Let us next describe these stability constraints. 



The most general Fierz transformed four-fcrmion interaction would include additional vector and axial-vector (qq) as well as scalar and 
pseudo-scalar, vector and axial-vector-like (qq)-interactions. However, in order to simplify the studying of the quark matter stability 
problem in the presence of external magnetic field, these additional terms are omitted. The role of the interaction in the diquark channel 
is discussed in the last Section 5. 
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Recall first that in the thermal and chemical equilibrium between vacuum and dense quark matter phases their 
(critical) chemical potentials, denoted as fj, c , are equal. Then, from the analogy with condensed matter physics, one 
can get the following relations between thermodynamical and dynamical parameters at which the two phases do 
coexist (called below stability constraints): 



mde 



(2) 



Here m vac is the dynamical mass of one-quark excitations of the vacuum and m^ense is the mass of quarks in the dense 
phase, i.e. inside the quark droplets (see also [T^l. 2 Note that the values of md ense , fj, c ,m va c are supplied by the 
thermodynamic potential f2(m; /i) of the system, where the average of the composite scalar cr-ficld m =< er > is the 
order parameter of spontaneous chiral symmetry breaking which in the thermal equilibrium is the global minimum 
point m a of the function f2(m; //) versus m. Clearly, m equals to the dynamical quark mass which is determined by 
the gap equation m — — 2G < qq >. (Obviously, if m = 0, the ground state is chirally invariant.) Hence, in order to 
find all the admissable values of parameters of the model, at which stable quark droplets could exist in the quantum 
system at moderate density, it is sufficient to solve the stability constraint with respect to these parameters. 
On this route one should perform several technical steps: i) get the thermodynamic potential f2(m;/j,); ii) find the 
global minimum point of f2(m; fi) vs m and establish its dependence on the external parameters of the theory thereby 
obtaining the phase structure of the model; iii) choose all the first order phase transition critical points fi c obeing 
the constraint J2J). 3 The corresponding values of allowed parameters obtained in this way are just the physically 
acceptable ones, at which stable quark matter might be described by the initial model. 
Now let us carry out this scheme on the basis of the chosen NJL model. 



B. Solution of the stability constraint at H = 



In the mean field approximation, the thermodynamic potential of the model (1) is given by the expression |13ll22ll24| : 



f2(?n; n) 



4G 



2N c N f 



d 3 p 



E. p + 6(n - Ep)(n - E p ) \, 



(3) 



where E p = \J m 2 + p 2 , N c is the number of colors, i.e. N c = 3, and Nf — 2. Using for the divergent integral in © 
a Lorentz noninvariant regularization (this procedure is simply the replacement of the infinite integration region in 
© in favour of the compact one: p 2 < A 2 , where the cutoff parameter A is an additional parameter of the theory to 
be fitted by the experimental data), one can find the following expression 



M(m;[i) = — 

N f 9(p - to) 
8^ 



3N f 



A(2A 2 + m 2 ) vro^+A 2 - m 4 In 



A + Vm 2 + A 2 



/z(2/z 2 - 5m 2 ) V V 2 - m 2 + 3m 4 In 



fj, + \/ H 2 — m 2 



■m 



(4) 



1. The vacuum case fi = 0. 



In this case, from Q one can get the following stationarity equation: 



d „ 3Nfm ( ir 2 . ■ 

^(rn;0) = ^[—-F(rn,A)}=«. 



where 



(5) 



F(m, A) = Ay/rn 2 + A 2 - m 2 In 



A + \Jm 2 + A 2 



(6) 



2 Indeed, one can interpret the chemical potential of any fermionic system as the lowest energy, needed for the particle in order to abandon 
the system. The equilibrium value fi c must be smaller than the mass of the vacuum one-particle exitations m vac (such particles have 
enough energy for leaving the vacuum, which is empty due to this reason). However, in order that quarks do not leave the dense phase, 
their mass "i^ense must be smaller than fi c - the minimum particle energy, necessary for leaving the system. 

3 At the critical line fi c of the 1st order phase transition, the order parameter m changes its value by a jump, i.e. at this line it is a 
discontinuous function vs input parameters of the model. At the critical line of the second order phase transition the order parameter 
is a continuous function, but first- or higher order derivatives of it are discontinuous. 
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The function F(m,A) monotonically decreases on the interval m S (0, oo) such that F(0,A) = A 2 and F(oo,A) = 
(F(m, A) = 2A 3 /(3m) + o(A/m) at m — ► oo). Hence, equation JSJ has a nontrivial solution m (G, A) = M ^ if 
and only if the constraint 3 ^ fG — F(0,A) < (or G > G cr n = 3] ^ fA 2 ) is fulfilled. In this case, the global minimum 
of the function fl(m; 0) lies at M. Thus, both the continuous and discrete chiral symmetries of the NJL-model are 
spontaneously broken. (At G < G cr it the stationarity equation Q has only the trivial solution m = 0. Thus, in this 
case the vacuum is chirally invariant, and fermions have zero mass. Since this phase is not adequate to the real QCD 
vacuum at T = p = 0, we will not consider in detail the undercritical values of the coupling constant in the present 
paper.) Evidently, among the three parameters A, G, M only two are independent ones at G > G cr it- In the following 
it is sometimes convenient to use A and M as independent free parameters of the model at G > G cr u . Clearly, M 
is here the mass of the one-quark excitation of the chirally-noninvariant vacuum with zero baryon density which, at 
H = 0, is just the quantity m vac from (J2J). Then, G = 3Nf p^ M A ) an d there is a one-to-one correspondence between 
G > G crlt and M > 0. 

2. The case p=fi0. 

Note first of all that for G < G cr n and arbitrary values of p, the chiral symmetry is not spontaneously broken down. 
However, above the critical coupling constant the situation is more complicated. Recall, that in this case the phase 
structure of the simplest NJL model with Nf = 1 and Lorentz invariant regularization was investigated in [25|- Using 
the same methods, it is easy to obtain the phase structure of model Q for the Lorentz-noninvariant regularization 
considered above. Thus, without any preliminary calculational discussions, we present the phase portrait of the NJL 
model (1) schematically on Fig. 1. On this figure one can see the (M, /i)-plane, which is divided into three domains: 
A, B, C. There is a symmetric dense phase with massless quarks for points (M, p) S A. In the region B, the phase with 
zero baryon density is realized. In this phase the chiral symmetry is spontaneously broken down, and quarks aquire 
a mass M. (In this sense, phase B models the real QCD vacuum.) Finally, the region C on this figure corresponds to 
the dense chirally noninvariant phase, in which quarks have a nonzero //-dependent mass. The boundaries between 
these phases are the critical lines of second order phase transitions (solid lines on Fig. 1), or first order ones (dashed 
curves on Fig. 1). In particular, the curve Ibc = : p = M} is the critical line of the second order phase 

transitions, since the second order derivative d 2 m / (dp) 2 of the global minimum point m of the potential fl(m; p) is 
a discontinuous function vs p on this curve. The curve Iab — {(M,p) : p = p c (M)} (all points of this critical curve 
lie below the line p = M) is defined by the equation 

fi(M; M ) = n(0;M), (7) 

i.e. on this line the order parameter of the system - the global minimum point m - changes its value by a jump from 
M to 0, and a first order phase transition takes place. Solving Q in the region p < M, one can obtain the manifest 
form of this curve: 

r q i i/4 

fi c (M) = \ 3AVA/ 2 + A 2 - 3A 4 - -M 2 F(M, A) \ , where M>M C = 0.56A. (8) 

Moreover, there are two tricritical points a and (3 on this phase portrait. Their coordinates are: a = (M,p), 
j3 = (M c , M c ), where M = 0.31A, p, = 0.37A and M c is defined in JHJ). 

Now one can easily see that constraints are fulfilled at M > M c (or in terms of coupling constants G and 
G crit , stability constraints (J2J are valid only at G > Gb ag = l.37G crit ). Indeed, for each such value of M there is 
a corresponding critical value of the chemical potential n c {M) © (recall, that < p c (M) < M) at which the two 
phases do coexist. In other words, in this case stable droplets of the dense phase A are surrounded by the chirally 
noninvariant vacuum (phase B) with zero baryon density. Obviously, the inequalities J5J are valid, since in the case 
under consideration rridense = and m vac = M . 5 The baryon density inside the quark droplet is easily calculated 
from the relation: 



Pdrop(M) = ' " { ' ' 



N c dp 



Nt 

(m = AvKM). (9) 



4 A point of the phase diagram is called a tricritical one if, in an arbitrary small vicinity of it, there are first- as well as second-order phase 
transitions. 

5 It might appear that points of the 1st order phase transitions critical line, which connects two tricritical points a and f3 of Fig. 1, are 
also suitable ones. However, in this case there is the coexistence of two dense phases, and the physical vacuum cannot be identified with 
one of them. So, we reject the corresponding values of M G (M, M c ) at H = 0. 
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Now, in the framework of the NJL model, one can imagine for M > M c the following quark matter transition in 
dependence on the average baryonic density p: At sufficiently small p a rather small amount of quark droplets, having 
a rather small size, appears in the space. The droplets are composed of massless quarks and surrounded by the 
vacuum with zero baryon density; however, inside each droplet the baryonic density equals to Pdrop(M) With 
growth of p the number of droplets as well as their sizes are increasing. When p reaches the value Pdrop(M), the total 
space is occupied by massless quarks. 

Finally, notice that the existence of stable quark droplets surrounded by a chirally noninvariant vacuum was also 
considered in the framework of a NJL model in 13]. However, there the problem was investigated only for three 
physically interesting sets of parameters: M = M X = 0.48A, A = 650 MeV (set 1), M = M 2 = 0.67A, A = 600 MeV 
(set 2) and M = M 3 = 0.88A, A = 570 MeV (set 3). Of course, the corresponding results are reproduced by our 
ones: for set 1 stable quark droplets are forbidden, but for sets 2,3 stable droplets are formed from massless quarks. 
In particular, it follows from J^J that pdrop — 2.86p for set 2 and pdrop — 4.19p for set 3 (here p is the ordinary 
nuclear matter density). In the present paper we have generalized these considerations and found the exact interval 
M > 0.56A = M c corresponding to stable quark droplet formation. 



III. PHASE STRUCTURE AT H 

Let us next study the phase structure of the NJL model JQ) containing u- and (i-quarks with electric charges 
e,\ = 2|e|/3 and e-i — — | e | /3 (e is the electric charge of electrons), respectively. Due to the difference in quark electric 
charges, the continuous chiral symmetry is violated in the presence of an external magnetic field, but the discrete one, 
q —>■ ij 5 q, is not. It is useful to consider first the vacuum case with p, = 0, where the thermodynamic potential of the 
system looks like: 

fi(m; 0, H) = g + p J™ J exp(-W) coth(| e ,iJ| s ), (10) 

where N c = 3. This expression can easily be transformed into the following one (see, e.g., |22j): 

Q{m;0,H) = n(m; 0) - £ (-1, Xi ) - \[x 2 - Xi ] lnx, + (11) 

j=i 

where 0(m; 0) is given in J2J at p = and Xi = m 2 / {2\eiH\) . Moreover, £'(— l,x)= d(,{v, x)/dv\ v —^x, where Q{v,x) is 
the generalized Riemann zeta- function |2(j. The thermodynamic potential 11011 - 111(1 is an even function with respect to 
to. This fact reflects the symmetry of the original model under discrete chiral transformation q — > i-f 5 q in the presence 
of an external magnetic field. Obviously, this symmetry might be spontaneously broken if the global minimum point 
to of f2(m; 0, H) is unequal to zero, but remains intact if m = 0. In the case G > G cr u, the global minimum point 
of 0(to; 0, H) is the solution of the stationarity equation (see p^ - ): 

JL —V(m; 0, H) = m{2F(M, A) - 2F(to, A) - I x (to) - J 2 (to)} = 0, (12) 
6 am 

where the function F(x, A) is defined in (JSJ, 

him) = \ ei H\{\nT( Xl ) - - ln{2n) + Xi - -(2 Xi ~ l)lnxj (13) 

and T(a;) is the first order Euler gamma- function [2(j. Below we suppose that H is a nonnegative quantity. 

The function F(x, A) vs x monotonically decreases on the interval x S (0, co) (see subsection II. B), and the functions 
Ii(x) vs x also monotonically decreases on the same interval from +oo to (see |22jV Hence, for arbitrary fixed values 
of H,M there exists only one nontrivial solution m (H) of equation (|12f) . which is the global minimum point of 
£l(m; 0, H). The quantity m (H) is a monotonically increasing function of H and m (H) — * M at H — > 0. It means, 
that at G > G cr n (or for all M > 0) the external magnetic field H enhances the discrete chiral symmetry breaking. 
Note that in the presence of an external magnetic field the parameter m vac from (j2Jl is equal to m (H), i.e. depends 
on the values of H. 6 



At G < G cr it, the magnetic catalysis effect of dynamical chiral symmetry breaking (DCSB) takes place. This means that at H = the 
NJL vacuum is chiral ly symm etric, but an arbitrary small value of the external magnetic field induces DCSB, and fermi ons acquire a 
nonzero mass m (H) |2 it l'27l (see also the reviews |2S| and recent papers on applications of the magnetic catalysis effect I29H V 
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Now let us consider the more general case, when H ^ 0,// 7^ 0. Using the methods of |22(, it is possible to obtain 
the thermodynamic potential: 



3\eAH 



0(m; fi, H) = fl(m; Q,H) - ^ ^ a k 9(fi - s ik ) < H\J - s} k - s-jt In 

z=l fe=0 I 



(14) 



where = ^/to 2 + 2|ei|-fffc, Qfc = 2 — Jot, and f2(m; 0, 7J) is given in (|llf) . The corresponding stationarity equation 
at G > G cr it looks like: 

^-^-n{m;fi,H) = m{2F(M,A)-2F{m,A)~h{m)-I 2 {m) + 
6 am 



^ |ei|gy^a fc g(^- s jfe )ln 



fe=0 



2 _ „2 



(15) 



It is useful to divide the (H, /x)-plane into regions Wfc, which have the following form: 



(H,n)={Ju k , where w fc = {{H, fi) : 2|ei|fl"fc < /i 2 < 2|ei|i?(fc + 1)}. 



(16) 



fe=0 



Then, the construction of the phase portrait is greatly simplified, since for each region ui n from lltifl only the terms 
with k = 0, .., n from the ei-sums in (|14[l - (|15[l are nonvanishing, and only the terms corresponding to k = 0, .., 2n + 1 
from the e2-sums in these expressions are nonzero as well. Hence, it is possible to perform step-by-step the numerical 
investigation of the phase structure of the NJL model in the regions too, u>\, etc. In particular, we have considered in 
detail the phase structure for the same three values of model parameter M = Mi, M2, M3 as in the paper [l3|. 

In the case M = Mi, we have studied the phase structure only in the region luq (116(1 . The phase portrait in terms 
of (H, /i) G uiq is presented schematically on Fig. 2, where L = { (H, /i) : /i = m (H)} denotes the upper estimate of 
the value fi c from @, i.e. m vac = m {H). Below L the (H, /i)-plane is divided into four domains A, B, C and Ci. 
The solid and dashed lines on this figure correspond to the critical curves of second and first order phase transitions, 
respectively. The points from the region A correspond to the massless, chirally invariant and dense phase. In region 
B there is the chirally noninvariant phase with zero baryon density. The one-quark excitations of the vacuum phase B 
have a mass equal to m {H), which depends on H and has no ^-dependency (the phase B can be identified with the 
real QCD vacuum without baryons). On Fig. 2 one can see also two another massive, i.e. chirally noninvariant, dense 
phases C and Ci as well as three tricritical points a,f3,j. Their coordinates are ^/2\e\Ha/A = 0.85, ^ Q /A = 0.45, 
y/2\e\Hi3/A = 0.82, [ip/A =0.47 and y/2\e\H^/A = 0.67, ^ 7 /A = 0.48. 

Now let us draw special attention to the boundaries Iab and Ibc of the vacuum phase B. The boundary Iab between 
phases A and B is the graph of the function fj, = [i c (H), where fi c (H) is the solution of the equation 



n(0;n,H) = Sl(m o (H);n,H), 



(17) 



i.e. on this curve the order parameter - the global minimum point of the thermodynamic potential - changes by a 
jump its value (from m (H) to 0), and a first order phase transition takes place. It turns out that the critical line 
Iab as a whole is arranged inside the a->o-region. However, the boundary Ibc-, which is also a critical line of first order 
phase transitions, where the order parameter changes by a jump from m (H) to some another definite nonzero value 
mi{H), is stretched through the regions ljq,.., wioo at least. Evidently, it tends to the point Mi (see Fig. 2), but 
to prove this hypothesis, we would need more calculational capacities. The curve Ibc is the graph of the function 
/i = fic(H), where p, c (H) is defined by the equation 



n{ mi (H); H) = Q(m (H);fi,H). 



(18) 



At M = M2,3 and H ^ the phase portraits at /i < m (H) is qualitatively the same as in the case H = 0: there 
are two phases A and B in this region, and the boundary Iab is the curve of first order phase transitions. 



IV. QUARK DROPLET STABILITY AT H 

In the previous section the phase structure of the NJL model was investigated at H ^ and for three particular 
values of the model parameter M: M = M 1 = 0.48A, M = M 2 = 0.67A and M = M 3 = 0.88A. Now, putting special 
attention just to these values of M , let us discuss the quark matter stability problem in the presence of an external 
magnetic field. Stability means that for fixed H there is the critical value /x c of the chemical potential at which the 
vacuum and dense quark phases coexist, i.e. at fi c a first order phase transition occurs. In this case, the value \x c 
must satisfy the stability relations 
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A. Magnetic catalysis of quark matter stability 

First of all, let us draw our attention to the case M = M\. The corresponding phase portrait is shown on Fig. 2, 
where points (H, \i = m {H)) are located on the curve L, and m (H) is the global minimum point of fi at fi = 0. 
Hence, m (H) is the mass of the one-quark exitations of the empty vacuum phase B and must be identified with m vac 
from © ■ The boundary of this phase is the first order phase transition line which consists of two curves Iab and I bc 
and lies below the line L. 7 So, at H > H a , there is coexistence of the zero density vacuum phase B with the dense 
phase A on the curve Iab (the one-quark exitations of the phase A have zero mass). Now, if one puts mdense = 0, then 
fi c = Hc(H) satisfies the stability constraints @, where fi c (H) is defined by the equation l|17l) . Hence, at H > H a , 
multiquark droplets which are schematically the bags of massless quarks, are allowed to exist in the vacuum. 

At H < H a and in the case (H, /i = jj, c (H)) € Ibc, where ft, c (H) is defined by the equation (|18J) . the function 
fi(m; /i, H) has two different global minimum points: the first one is located at m = m (H), the second one at some 
definite nonzero point m\{H). So, at H < H a there is a coexistence of the vacuum phase B with the dense phase 
C on the curve Ibc (the value mi(H) is just the mass of the one-quark exitations of the ground state of the phase 
C). Now, if one puts mdense = m i(H), m vac = m (H), then \i c = fi c (H) satisfies the stability constraints (JSJ. This 
conclusion can be obtained from a detailed analysis of the stationarity equation l|15f) . At A = 650 MeV it is also 
confirmed by numerical calculations depicted on Fig. 3 showing graphs of the functions m vac , fi c and mdense vs H 
corresponding at y / 2\e\H/A < 0.85 to those points of the line Ibc, which lie inside the region ojq. So, at H < H a 
stable quark droplets composed of massive quarks are allowed to exist in the NJL model. 

Recall, that in the case M = Mi the stability of quark droplets is absent at H = (see Section II as well as [13J). 
However, as we have proved, in the presence of an external magnetic field there might exist stable quark droplets, i.e. 
an external magnetic field promotes the stabilization of quark matter. Moreover, it follows from our consideration 
that at H ^ 0, in contrast to the case H = 0, stable multiquark clumps might be composed both from massless quarks 
(at H > H a ), and massive ones (at H < H a ). These are the main results of our paper. 

Remark, that usually the value of the cutoff parameter A varies in the interval (0.5 -j- 1) GeV. In this case, the 
magnitude of H a can be estimated to be of the order 10 19 Gauss. 8 According to modern estimates, the strength of 
the magnetic field at the surface of neutron stars can reach 10 15 Gauss, and at the star core it may be as large as 
10 18 Gauss [s^. Hence, for M — Mi as well as for some neighbouring values of M, the NJL model may predict the 
existence of stable multiquark clumps, composed of massive quarks, inside neutron stars. 

Note finally, that for arbitrary values of M £ (M,M C ) (see Fig. 1), quark droplet stabilization is induced by the 
external magnetic field in the same manner, as in the case M = M\. 

At M — M2.3 and in the presence of an external magnetic field, the zero density vacuum phase B (rn vac = m (H)) 
might coexist only with the dense phase A (mdense = 0) at \x — /i c . Critical values \x c of the chemical potential in 
these cases are depicted on Figs 4, 5 for M = M2 and M = M3, respectively. It follows from these figures that the 
stability constraints © are also fulfilled, so the stable quark droplets are allowed to exist. As in the case H = (see 
at M = Af2,3 and H 7^ the finite lumps of stable quark matter are schematically the bags of massless quarks 

only. 

The phase structure of the NJL model at H ^ and fi ^ was studied in 0, ^ , but for qualitatively other 
values of model parameters. In particular, the stability problem of multiquark clumps was not considered in these 
articles. Thus, in view of our present results, let us now discuss it in these cases, too. In |2l|. in connection with 
the magnetic catalysis effect, the influence of an external magnetic field on the dense NJL model was investigated at 
G < G cr it- It turns out that at sufficiently small values of /z and arbitrary values of H the external magnetic field 
induces the rearrangement of the model vacuum: at H = it is chirally invariant, but the chiral symmetry breaking 
is induced at H ^ (magnetic catalysis effect [2jJ,|28j). Due to the zero baryonic density it is the vacuum phase B. 
At larger values of \x the dense chirally symmetric phase A with massless quarks is realized. As shown in [21| , a first 
order phase transition occurs between the A and B phases. Now one can easily check that the stability condition @ 
of the present paper is fulfilled at G < G cr u, H ^ 0, too. Hence, in this case stable droplets, composed from massless 
quarks, are allowed to exist. These droplets are surrounded by the chirally noninvariant vacuum. In papers |2^| the 
parameter M was fixed in such a way, that M < M < M c (see Fig. 1 of the present paper), i.e. at H = the quark 
matter is not stable in this case. Then, a rather complicated phase portrait was shown to exist in the (H, /i)-plane. 
A more detailed analysis results in the statement that the stability constraints @ can be solved for these values of 
model parameters as well (in a similar way, as in the case M = Mi). 

Hence, at H 7^ 0, in contrast to the H = case, stable quark droplets are allowed to exist in the framework of NJL 
model (1) for all values of the coupling constant G > 0: at G < Gbag this property is induced by an external magnetic 
field, whereas at G > Gbag it is an own property of the model. 



7 On Fig. 2 there is another critical line, where 1st order phase transitions take place as well. It is the line which connects tricritical points 
fi and 7. However, it is the boundary between two dense phases, and no one of them might be identified with the physical vacuum. 

8 This estimate can be obtained from the known relation m^,/\e\ = 4.4 X 10 13 Gauss (where m e and e are the mass and the electric charge 
of an electron, respectively), which permits to connect (MeV) 2 and Gauss units. 
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B. Energetical approach to the stability problem 

It might seem that at M = M2,3 and H ^ the solution of the quark matter stability problem is qualitatively the 
same as at H = 0. However, there are some peculiarities which can be found by an energetical consideration of this 
problem. Analogously to the case at H = [llj, let us define the energy density e at H ^ by the following relation 

e = n(m; n, H) - 0(m o (if); 0, H) + N c pp, (19) 

where p = —dQ/(N c dp) is the baryonic density and m (H) - the dynamical quark mass in the vacuum phase B. A 
very important physical quantity is e/p - the energy per one baryon. At H = the graph of this function vs p is 
presented on Fig. 6 for M = M% and A = 600 MeV (see also E3)- The absolute minimum of this function lies at the 
point p — pdropiMz) = 2.86po which is equal to the baryonic density inside stable quark droplets. At the point p = 
there is a local minimum of e/p, corresponding to the real empty physical QCD-vacuum, i.e. to the phase B. Due to 
the configuration of the curve e/p vs p at H = 0, one can conclude that the vacuum, i.e. the state with p = 0, is 
a metastable one. It means that not for all values of baryon density fluctuations the system will pass automatically 
to the mixed phase with stable quark droplet formation. Indeed, the curve e/p has a local maximum at some point 
Pmax- So, if a quark droplet with p < p max is created in the vacuum, it is energetically favorable for the system to 
reach zero density, i.e. to let tend the volume of the quark droplet to infinity thus pushing out all quarks from the 
system. As a result, one gets the initial empty vacuum. However, if a quark droplet with p > p max is formed, then it 
will be a stable one, since in this case the droplet simply changes its volume in order to reach the most energetically 
favorable density PdropiMi). 

The situation is essentially different in the presence of an external magnetic field. In this case the graph of the 
function e/p has a local maximum in the point p = (see Fig. 6), so at arbitrary small baryonic density fluctuations 
the stable quark droplets will appear in the system. A similar property of the function e/p is observed for the case 
M = Mi (see Fig. 7) as well as for other values of the parameter M. Hence, the external magnetic field promotes 
the creation of the stable quark matter. 

Let us denote the global minimum point of the function e/p as pdrop (this value is the baryonic density inside the 
stable quark droplets). Then, the binding energy per quark Eb inside stable droplets is given by the relation 

- m vac = m vac - p c , (20) 

P=Pdro V ) 

where m vac = m (H). Of course, both quantities Eb and pdrop depend on H. Their graphs are presented on Fig. [S] 
and Fig. respectively. Recall, that in the case M = Mi the quark binding energy Eb is about 20 MeV at H = 
(see E3)- Now one can easily see from Fig. [S]that at nonzero H this quantity might be several times larger, than at 
H = 0. So, the external magnetic field enhances the stability of quark matter. 

Note that on Fig. [5] the density curve corresponding to the case M = M^ has a cusp at the point \j2\e\H /Awl. 
The matter is that in the chirally symmetric phase such quantities, as particle density, magnetization etc, have a 
discontinuity of their derivatives over p and H, when a new Landau level is switched on [22j . In the case under 
consideratio n in a ddition to the lowest Landau level the first one for the d-quarks is brought into account just at 
the point y/2\e\H/A w 1, so the cusp for the density curve appears. (On Fig. 03 a cusp on the curve Eb vs H also 
appears at the point y/2\e\H / A « 1 in the case M = M2, which is, however not so sharp. One can easily see a cusp 
on the curve p c vs H of Fig. 4 as well.) Of course, the smaller the lowest value of the H- interval which is taken 
into account, the more cusps will appear on the density curve as well as on the graphs of other physical quantities. 
At H — ► there are infinitely many cusps on these curves, which results in the experimentally observed oscillating 
phenomena of corresponding physical quantities (see p^ - ). 

V. SUMMARY AND DISCUSSIONS 

In the present paper the role of an external magnetic field H on the formation of stable quark matter droplets inside 
the chirally noninvariant vacuum is investigated in the framework of the simple N JL model . Earlier, the problem 
of quark matter stability was considered in terms of more general effective theories with four-fermionic interactions 
as well, but only at H = [H El EH ■ 

Thus, we began our study with the case of a vanishing extern al m agnetic 
field, too. Moreover, we used the same thermodynamic approach for treating multiquark droplets as in |13|. In this 
case, any multiquark clump is imagined as a droplet of the dense quark phase of the system, which is surrounded 
by a chirally noninvariant vacuum phase, and there is thermal and chemical equilibrium between both phases. This 
qualitative picture can be transformed into a quantitative constraints on dynamical quark masses and the chemical 
potential at which the coexistence of two phases is possible. Moreover, it was shown in EH that a t some values of 
the physically acceptable model parameters (coupling constants, cutoff parameter, etc) the stability constraints (0) 
are not fulfilled for the NJL model at H = 0. In those cases stable multiquark droplets are forbidden, and the NJL 
model is not applicable for the description of dense quark matter. 



E B = - 
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In our paper we have in the first step obtained the general solution of the stability constraints J2Jl at H = for 
the NJL model considered above. Here, we have considered only the physically interesting case, when G > G cr n — 
ir 2 /(3NfA 2 ), in which case chiral symmetry is spontaneously broken down and quarks acquire a dynamical mass 
M (at jU = 0). Due to the stationarity equation (JSJ, the quantities M and G (at G > G cr u) are in a one-to- 
one correspondence. Thus, we find it convenient to use the value M (instead of G) as a model parameter (if M 
changes from to oo, then G monotonically increases from G cr u to oo). Stable quark matter droplets inside the 
chirally noninvariant vacuum are then allowed in the model under consideration only at M > M c = 0.56A, i.e. at 
G > Gbag = l-37G cr ii, and they are composed of massless quarks. Notice also that in 13] the stability problem 
was investigated only at three physically admissable sets of parameters: M = Mi = 0.48A, A = 650 MeV (set 1), 
M = M 2 = 0.67A, A = 600 MeV (set 2) and M = M 3 = 0.88A, A = 570 MeV (set 3). There it was pointed out that 
for set 1 stable quark droplets are forbidden, but for sets 2 and 3 stable quark matter is allowed to exist. Of course, 
these results coincide with our ones at H = 0. 

Next, we took into account an external magnetic field and studied the problem of quark matter stability in detail 
at the same values of the parameter M as in ref.0]. It turns out that at any value of H in the cases M = M 2 ,3 stable 
quark droplets are allowed to exist. Moreover, as at H — 0, these multiquark droplets consist of massless quarks. 9 
The most interesting results were obtained in the case M = M\ , where we have found that an external magnetic field 
induces stable quark matter formation. Indeed, at H = multiquark droplets are not stable in this case. However, 
at H they are allowed to exist. Furthermore, there exists a value of magnetic fields, H a ~ 10 19 Gauss, such that 
at H < H a quark droplets are formed by massive quarks. At H > H a , multiquark droplets are composed of massless 
quarks. We have also checked that in a similar way, for other values of M from the interval M £ (M,M C ) (see Fig. 
1), stable quark droplet formation is also induced by the external magnetic field. 

Note that one can even assert that at H ^ stable quark droplets are allowed to exist in the initial NJL model at 
arbitrary values of the coupling constant G > 0. At G > Gbag it is the native property of this theory even at H = 0, 
but at G < Gbag the stability of the quark matter is induced (catalized) by the external magnetic field. 

Moreover, using an energetical approach, we have shown that an external magnetic field promotes the creation as 
well as enhances the stability of quark droplets at G > Gbag- 

Finally, let us add some remarks about the role of color superconductivity (CSC) in the quark matter stability 
problem considered above. This phenomenon was predicted quite recently in the framework of instanton-induced- as 
well as NJL-type models with four - fermionic interactions (see the reviews 31] or some of the recent papers 32, 33] 
and references therein). At zero temperature it is proved that CSC must exist at asymptotically large values of the 
chemical potentialfor baryonic densities), however, the lower boundary for this effect is not established up to now. 
In particular, in [34J CSC was investigated just in a slightly extended version of the NJL model (1) at neighboring to 
the set 1 values of parameters, i.e. for M « Mi and A = 650 MeV, but with an additional four-quark term included, 
which is responsible for the interactions in the diquark channel (with a corresponding coupling constant G ). It was 
realized at H = 0, that usually with growing value of /x, the chiral invariance of the system is restored before or just 
at the CSC phase transition. Moreover, at fixed G, the parameters of the CSC phase transition strongly depend on 
the value of n = G /G, while the features of the chiral phase transition do not seriously depend on n. For example, in 
the NJL model (1) at small values of k < 1/4 the transition to the CSC phase occured at [i ~ 700 MeV, but at K ~ 1 
it is at n ~ 300 MeV, whereas for wide interval of K-values (G is fixed) the chiral symmetry is restored at fi ~ 300 
MeV H|. 

For simplicity, we have put G = in our consideration of the quark matter stability problem at H 0. In the 
light of the above discussion we believe that our results are not qualitatively changed at least for small values of k, 
i.e. in the case, where CSC occurs for rather large values of /i ~ 700 MeV [34j • In fact, in this case CSC can not 
seriously influence stable quark droplet formation which takes place at values of the chemical potential much less 
than 700 MeV (see the graphs of \x c on Figs 3-5). At present time there are several indications on the possibility for 
CSC to take place only at chemical potential values which indeed are outside the scope of our consideration, i.e. at 
/i >> 300 MeV: i) the CSC phenomenon was not discovered up to now in heavy ion collision experiments, ii) It was 
argued in |35| that two-flavor CSC is absent inside quark cores of neutron stars, where baryonic density is several 
times the density p a of ordinary nuclear matter, iii) In j3|| it was shown that due to the color Meissner effect as well 
as due to the presence of background chromomagnetic fields (gluon condensate) the CSC is shifted to higher densities 
(chemical potentials), iv) As in the case with ordinary superconductivity, the sufficiently strong external magnetic 
field prevents a system from CSC as well etc. In these cases the ordinary, i.e. without CSC, quark matter might be 
formed at rather low densities. To describe such a processes one could use models of the type (1) with zero or very 
small value of the diquark channel coupling constant G . 

Recently, we have studied already the magnetic catalysis of the quark droplet stabilization in the scope of a simpler 
one-flavor NJL model [3^ . Next, we are planning to consider the quark matter stability problem in more realistic 



9 It is necessary to note, that for such values of M which are greater, but quite near to M c JSJ (for example, at M c < M < 1.01M C ), 
there exists a finite interval of the external magnetic field values, where stable quark droplets are formed by massive quarks. But for 
values of the external magnetic field outside of this interval, stable multiquark droplets are composed of massless quarks. 
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NJL models including vector interactions as well as the third flavor of quarks etc. 
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Figure Captions 



Fig. 1. Schematic phase portrait of the NJL model at nonzero /x and for arbitrary values of the quark mass M 
at \i = (the case G > G cr it). Phases B and C are massive nonsymmetric phases, A is a chirally symmetric and 
massless phase. In the phase B the quark mass is equal to M. In the phase B the particle density in the ground 
state is equal to zero, whereas in phases A, C it is nonzero. Solid and dashed lines represent critical curves of second 
and first order phase transitions, respectively; a and (5 denote tricritical points. In particular, the curve Ibc has the 
following form: Ibc = {(M /i) : \i = M}, the curve Iab is defined by the equation (JTJ), M — 0.31A, M c = 0.56A. 
Only the values /i c (M) (see of the chemical potential, corresponding to points of the critical line Iab, obey the 
stability constraints J2J). 

Fig. 2. Schematic (ff, ^)-phase portrait of the NJL model at M = M x = 0.48A (M < M x < M c ), where A is 
the chirally symmetric massless phase, B, C and Ci are chirally noninvariant massive phases. Ground state baryonic 
density is zero in the phase B, and nonzero in phases A, C, Ci. Dashed and solid lines are critical curves of the first- and 
second order phase transitions, respectively. (Iab is defined by the equation ^17$ ) The curve L={(H, fx) : fi = m (H)}, 
where m (H) = m vac is the quark mass at fx = 0. fi\ = 326 MeV is such a value of the chemical potential that the 

point (Mi,/ii) lies on the curve a/3 of Fig. 1. The external magnetic field values, corresponding to the tricritical 
points a,/3,7 are the following: y/2\e\H a /A = 0.85, y / 2\e\H fj /A = 0.82 and y / 2\e\H~ f /A = 0.67. Only the points of 
the critical lines Iab and Ibc obey the stability constraints lf2"|. 

Fig. 3. The behaviour of quantities m vac = m (H) 1 /i c and mdense vs H which characterize the coexistence of 
the phase B with A and C ones in the case M = Mi and A = 650 MeV. If ^ l /2\e\H /A > 0.85 these parameters 
describe the first order phase transition (curve Iab) between A and B (in this case mdense = and \x c = fi c (H), where 
/x c (H) is the solution of the equation fT7| ). if ^2\e\H/ K < 0.85 - the phase transition between B and C (in this case 
mdense = TOi(£f) and fi c = jl c (H), where jl c (H) is defined by the equation l|18|l 'l. Only the part of Ibc, lying in the 
region ujq is considered. These values of m vac , \i c and mdense are the solutions of the stability constraints 

Fig. 4. The same as on Fig. 3, but for M = M 2 and A = 600 MeV (in this case mdense = 0). The dotted 
horizontal line corresponds to the critical chemical potential value /i c at zero external magnetic field. For these values 
of m vac , \i c and mdense the stability constraints @ are fulfilled. 

Fig. 5. The same as on Fig. 3, but for M = M3 and A = 570 MeV (in this case mdense = as well). The dotted 
horizontal line corresponds to the critical chemical potential value /i c at zero external magnetic field. For these values 
of m vac , /i c and mdense the stability constraints (J2J) are fulfilled. 

Fig. 6. The energy per baryon in the NJL model (1) at M = M 2 , A = 600 MeV for H = and y/2\e\H/A = 1.04. 
The solid line corresponds to the massive solutions of the stationarity equation, the pointed line to the massless one. 



Fig. 7. The energy per baryon in the NJL model (1) at M = Mi, A = 650 MeV for v / 2\e\H/A = 0.82 and 
y/2\e\H / A — 1.04. The solid line corresponds to the massive solutions of the stationarity equation, the pointed line 
to the massless one. 

Fig. [8| The binding energy per quark for a stable droplet in the NJL model (1) vs external magnetic field in two 
cases: M = Mi, A = 650 MeV (solid line) and M = M 2 , A = 600 MeV (pointed line). 

Fig. |9| Baryon density pdrop inside a stable quark droplet in the NJL model (1) as a function of external magnetic 
field. The solid line corresponds to the case M = Mi, A = 650 MeV, the pointed one to the case M = M 2 , A = 600 
MeV. 
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